Kinks in the Hartree approximation 
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The topological defects of A<^'* theory, the kink and antikink, are studied in the Hartree approx- 
imation. This allows us to discuss quantum effects on defects in both stationary and dynamical 
systems. The kink mass is calculated for a number of parameters, and compared to classical, one 
loop and Monte Carlo results known from the literature. We discuss the thermalization of the sys- 
tem after a kink-antikink collision. A classical result, the existence of a critical speed, is rederived 
and shown for the first time in the quantum theory. We also use kink-antikink collisions as a very 
simple toy model for heavy ion collisions and discuss the differences and similarities, for example in 
the pressure. Finally, using the Hartree Ensemble Approximation allows us to study kink-antikink 
nucleation starting from a thermal (Bose-Einstein) distribution. On a qualitative level, our results 
show only few dissimilarities with the classical results, but on a quantitative level there are some 
important differences. 
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I. INTRODUCTION 



In recent years the increasing numerical power of com- 
puters has made it possible to follow quantum fields in 
real time, studying phenomena such as thermalization, 
phase transitions and the formation of topological de- 
fects. Until recently this was mainly done using the clas- 
sical approximation, which can be justified when occupa- 
tion numbers are large. However, this is never the case 
for all momentum modes, and for systems far from equi- 
librium, this approximation may also break down. The 
Hartree and large-n approximations do capture quantum 
effects but, until recently, the available numerical power 
did not allow for the simulation of inhomogeneous sys- 
tems, thereby excluding the study of many physically 
interesting phenomena such as topological defects. It 
also has problems with describing thermalization, due to 
the lack of direct scattering. Recently, approximations 
based on an expansion of the 2PI effective action (see 
Ref. U and references therein) and also on a truncation 
of Schwinger-Dyson equations (see Ref. f5| and references 
therein), which do include direct scattering have been 
shown to lead to a quantum thermal equilibrium. Unfor- 
tunately, these methods are still numerically expensive 
and are not yet able to describe inhomogeneous systems. 
In two papers 0, we introduced the Hartree ensem- 
ble approximation, which tries to combine the good fea- 
tures of both the classical approximation, which is able 
to thermalize (classically), and the Hartree approxima- 
tion, being a quantum approximation and thus free from 
the Rayleigh-Jeans divergencies. In the Hartree ensem- 
ble approximation, the initial density matrix is expanded 
on a (overcomplete) basis of Gaussian density matrices. 
Expectation values can then be written as a weighted 
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average of Gaussian expectation values. The equations 
of motion are derived using the Hartree approximation 
in each of the "realizations." The partial "mean fields" 
of these different "realizations," i.e., the expectation val- 
ues of the field in the different Gaussian initial states, are 
typically inhomogeneous and the interaction of the quan- 
tum modes with the Fourier modes of these "mean fields" 
leads to approximately thermal distributions. However, 
since these "mean fields" are typically inhomogeneous, 
topological defects can and will be present if the symme- 
tries of the theory allow for them. These will be the topic 
of this paper, we will look at the (j)'^ theory in 1-1-1 di- 
mensions and its defects, kinks, and make a comparison 
between the classical theory and the Hartree approxima- 
tion. 



The rest of the paper in organised as follows. We will, 
in the remainder of this introduction, briefly review the 
equations of motion and our definition of instantaneous 
particle number, used in studying the distribution func- 
tion. In Sec. ^ we will derive the classical kink solutions, 
as this will facilitate comparison with the Hartree ap- 
proximation. Then, in Sec. IIIII we will discuss Hartree 
kink-antikink pairs, initially at rest and compare with 
the classical approximation. A study of the Hartree kink 
mass will be made. In Sec. IIVI kink-antikink collisions 
will be studied. We will look at the critical speed: be- 
low a certain speed the pair annihilates, while above they 
bounce back. This has been extensively discussed for the 
classical theory, but not yet for the quantum theory. We 
will study whether the system thermalizes after the col- 
lision and use the collisions as a toy model for heavy ion 
collisions, looking at the central region and its depen- 
dency on coupling and initial speed. As a final topic we 
will, in Sec. look at thermal kink nucleation, starting 
from a Bose-Einstein thermal ensemble, as in Ref. Q. 
The conclusions are in Sec. IVII 
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A. Equations of motion 



B. Observables: Particle number 



In this section we very briefly discuss the Hartree equa- 
tions of motion for the (j)'^ theory, as used in the separate 
"reahzations" of the initial ensemble. For more informa- 
tion see Refs. [11113. 

The Lagrangian density for our 1+1 dimensional sys- 
tem is given by^ 



(1) 



2" " 4 

giving rise to the Heisenberg (operator) equations of mo- 
tion 



(2) 



By introducing expectation values of products of field op- 
erators one obtains a hierarchy of equations: the equation 
for the n-point function will contain also the (n -I- 2)- 
function. Unless the initial density matrix, which is used 
in the calculation of the expectation values, is Gaussian, 
this hierarchy will be infinite. The Hartree approxima- 
tion truncates it self-consistently by imposing Gaussian- 
ity, thus factorizing all n-point functions with n > 2 into 

1- and 2-point functions. It is convenient to expand the 

2- point function in a set of mode functions, which in the 
free theory are simply plane waves. This expansion is 
most easily implemented in the field operator itself^ 



{X, t) = 4>{X, t) + ^ Ja{x, t)ba + h.C. 



(3) 



It is important to note that, due to the Hartree approx- 
imation, the creation and annihilation operators and 
6jj can be chosen time- independent. They also satisfy 
the usual commutation relations. The complex mode 
functions fa{x,t) form a complete and orthogonal ba- 
sis with a Klein- Gordon- type inner product. However, 
the a-label in general does not coincide with the Fourier 
label in time-dependent systems. 

The Heisenberg equations of motion in the Hartree ap- 
proximation can thus be written as 

4> = A (j) - [n'^ + \(jy^ + 3AC]0, (4a) 

h = A/a - + + 3Aq/„, (4b) 



with 



C^Y.^2nl + l)\U 



(4c) 



In 1-f 1-dimensions, only the mass needs renormalization, 
due to the divergence of the mode sum C{x,t). We will 
do this by absorbing its divergent vacuum part in /i^. In 
terms of the renormalized /z^^j^ we have obtained a finite 
theory. In the rest of the paper we will only consider 
states for which n9, = 0. 



^ We use the metric (—1, 1). 

^ We assume a large finite lattice size L and use finite size notation, 
i.e., a discrete spectrum. 



In order to study thermalization we use the same def- 
inition for particle number as before ^ ^ 13 1 using the 
2-point functions. 



S{t,x - y) ^ {<P{t,x)^{t,y)} ~ {cf>{t,x)) m,y)), (5a) 



U{t,x - y) = {f{t,x)T:{t,y)) - {iT{t,x)) (7r(i,y)), (5b) 

where the overline denotes an average over the center 
of mass coordinate and optionally a time interval and/or 
ensemble of initial conditions: 



{<P{t,y + z)(t){t,y)) = 
N 2^ T„ 



AT^TXI / dym',y + z)cf>(t',y)), (6) 

^^.^^ Jt-s/2 Jo 



where N is the number of ensemble members, L the lat- 
tice size and S a small time interval. In principle the 
(asymmetric part of the) 07r correlation function could 
also be used, but as it vanishes in equilibrium, we will 
refrain from doing that. 

We can split the 2-point functions in the separate con- 
tributions from the mean field and the mode functions: 



S{t,x - y) ^ S'{t,x - y) + S'^{t,x - y), 



(7) 



where 



S'{t, x~y)^ <j>{t, x)(b{t, y) ~ cj,{t, x) cj,{t, y), (8a) 
S^it,x-y) = Cit,x;t,y), (8b) 

and similarly for U{t, x ~ y). These expressions only de- 
pend on the relative coordinate x—y. We take the Fourier 
transform 



Sk(t) 



1 



dz e 



'ikz 



Sit,z), 



(9) 



with k = (0, ±1, ±2, • • • )(27r/L). Analogously to their 
expressions for a free theory we can use them to define 
instantaneous particle numbers: 



Sk{t) 
Uk{t) 



nk[t) + - 



1 



Ukit) 



LUk{t) ' 
UJkit). 



(10a) 
(10b) 



From these relations we can extract nfe(t) and cokit) as 
a function of time. Compared to other definitions of par- 
ticle number, this one has the advantage of being man- 
ifestly real, thus making them particularly suitable for 
studying out-of-equilibrium systems. 
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II. CLASSICAL KINK SOLUTIONS 

When looking for stable solutions to the classical equa- 
tions of motion of the (p'^ theory in an infinite volume, the 
field has to be in the vacuum at spatial infinity in order 
to keep the energy finite. In the "broken phase" this still 
leaves some important freedom. For the trivial solution 
both ends will be in the same vacuum. However, there 
is also a nontrivial solution, where the field is in the dif- 
ferent vacua in the different spatial infinities and hops 
over the central barrier somewhere in between. This is 
the so-called (anti)kink solution. The exact location of 
the kink does not influence the energy and an infinite set 
of degenerate solutions exists, but transforming the kink 
solution into the trivial solution would cost an infinite 
amount of energy, and hence the topology has rendered 
it stable. See Ref. for a standard introduction into 
the topic. 

It is possible to calculate this lowest non-trivial solu- 
tion using a very elegant and general method, the Bogo- 
mol'nyi equation In the broken phase the classical 
static Hamiltonian (energy) is given by 



E = 



— dx ■ 



A.l 



(-03 _ y^^) 



2 ^3 



(11) 



-oo) 



where (/)' = dx<t>- The boimdary term is equal to a finite 
constant and determines the topological sector, trivial or 
nontrivial. For given boundaries the energy is therefore 
bounded from below by the first term which is larger 
than, or equal to zero, the so-called Bogomol'nyi bound. 
The lowest energy state can be found by solving the sim- 
ple first order differential equation 



which has solutions 



±w tanh ^-j{x — xq) , 
±v. 



(12) 



(13) 



Here xq is the point where 4> goes through and m = 
v2At^ is the quasiparticle or "meson" mass in the broken 
vacuum. In the top line, the -I- solution corresponds to 
a kink, the — solution to an antikink. The width of the 
kink is inversely proportional to the mass. 

The kink energy density can be found from the top line 
of using ^ 



^kink 



2\2 



1 



■A 



Aw'* tanh^ 



'-{x - xo) 



1 



(14) 
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FIG. 1: Classical kink solution Hl.'ifl and its energy density 
11411 as a function of x in units . 



Both the solution H13|l and the energy density H14|l are 
plotted in Fig. ^ 

The kink mass, defined as the energy of the static kink 
or antikink, can be found either by integrating II14I) or 
directly from the boundary term in Hll(l 



kink 



m 
3A 



(15) 



Note that it diverges in the limit A ^ 0, showing the 
nonperturbative nature of the configuration. 

From the static kink solution H13(l we can easily find a 
moving kink solution by a Lorentz transformation: 



4>{x, t) = V tanh 7-^(2; 



Xq — ut) 



(16) 



where 7 = l/Vl — u^- Note that this is not a solution of 
the Bogomol'nyi equation H12(l . since it is not static, how- 
ever it is a solution to the manifestly Lorentz covariant 
field equations. The total energy of this moving solution 
is equal to -Ekink = 7-^kink, as expected for a relativistic 
particle. In summary, a Lorentz boost has the effect of 
increasing the kink mass while decreasing its width. 



A. Lattice (arte)facts 

We end this section by mentioning some consequences 
of the use of a finite lattice. 

The topological sector is determined by the type of 
boundary conditions: when imposing periodic boundary 
conditions the net kink number is 0, but a kink-antikink 
pair can be considered. Although such a configuration is 
inherently unstable, its survival time depends heavily on 
the relative distance, as a result of the exponential at- 
traction between them. When using antiperiodic bound- 
ary conditions, the kink number is ±1. By shifting the 
kink over 1 lattice size L through the boundary, the kink 
number changes sign. This means one cannot distinguish 
between kink and antikink. However, the defect itself is 
stable, as expressed by the nonzero kink number. 

An unwanted side effect of the use of a lattice is radi- 
ation, resulting from discretization errors and interfering 
with the kinks themselves. Some solutions have been 
proposed, see for example Gleiser and Sornborger p] and 
Speight and Ward H El [ij- In Ref. Q a damping 
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FIG. 2; Kink-antikink initialisation procedure. 



method is used to remove the radiation before it can in- 
terfere. Speight and Ward H ES 113 propose a lattice 
discretization preserving the Bogomol'nyi bound, thereby 
also removing the radiation result ing from the discretiza- 
tion. However, according to Ref. [iJl this latter solution 
is only suitable for free kinks, while in dynamical systems 
a standard discretization would be simpler with effec- 
tively the same accuracy. We will not bother too much 
about this issue. If it will be necessary to accurately 
follow an approaching pair over a long period of time, 
large lattices will be used and the kink and antikink will 
be put close together initially: the radiation is mainly 
emitted backwards and therefore reenters via the peri- 
odic boundary condition. Thus, by setting the soliton 
pair close together, we can let it take a long time for the 
radiation to reach them again. 

Finally, the use of (anti)periodic boundary conditions 
in a finite volume introduces a small mismatch in the 
derivative of the field at the boundaries, especially if the 
kink and antikink are close to them. We therefore use 
the initialization setup sketched in Fig. [2 we take sev- 
eral configurations in several adjacent lattices and use 
the resulting configuration in one of them. 



'I^KKi^^'^)] = -V+ 

I t—0,aQ 



f 

V < tanh 

^ — _ ^ 



tanh 



^—[x+[xo — ut\+n- L) 



III. p , , 

7 — [x — [xo — ut\ + n ■ L) 



(17) 



A. Initial condition 

Unlike in the classical theory, we cannot use the Bogo- 
mol'nyi argument to derive a static soliton solution to the 
Hartree equations However, if the coupling is not too 
large, or more precisely, if the two vacua are well sepa- 
rated and the barrier between them is high, the quantum 
corrections will be relatively small. Furthermore, away 
from the physical position of the kink, the field resides in 
one of the vacua, where an exact solution of the Hartree 
equations is known. Therefore using the classical kink so- 
lution as initial condition for the mean field, while using 
the free field plane wave solutions for the mode func- 
tions will be close to a Hartree kink solution, provided 
the coupling is not too close to the phase transition, i.e., 
the vacuum expectation value v is (substantially) larger 
than 1. After the configuration has evolved, using the 
Hartree equations of motion, from such an initial con- 
dition, it will, at least for a while, oscillate around the 
actual stationary solution. 

A way to improve on this initial condition is to add a 
damping term 



(18) 



to the right hand side of the mean field equation^ (|4a|) and 
evolve from the aforementioned initial condition to obtain 
an approximately stationary solution, which can then be 
used as a new initial condition. We will discuss this setup 
further in Sec. IIVI when considering colliding kinks. In 
Sec. IIII CI we will use a damping term to determine the 
quantum kink mass. 

Another possibility would be to start with the initial 
conditions as proposed in Ref. However, these still 
are not stationary under the equations of motion and the 
advantage is therefore only minor, so we choose to use the 
simpler vacuum mode functions. 

The numerical results presented in the next sections 
are all obtained using the simple initial condition, with 
mode functions of free field form and Eq. (|17(l for the 
mean field. 



where N is typically around 8. 



B. Numerical results: Static kink decay 



III. HARTREE KINKS AT REST 

Although interesting in itself the classical theory is not 
the full story and it is therefore important to study quan- 
tized kink solutions. In order to study dynamical pro- 
cesses, such as collisions, one has to use approximation 
schemes, which are nonperturbative, real-time and which 
can handle inhomogeneous configurations. For this pur- 
pose the inhomogeneous Hartree approximation as stud- 
ied in Ref. and Q should be particularly appropri- 
ate, being a nonperturbative semiclassical approximation 
scheme. 



In this subsection we are discussing the evolution of 
a kink-antikink configuration, initially at rest at a max- 
imum distance in the periodic volume. A comparison 
will be made between Hartree and classical dynamics, 
using identical initial conditions, in the sense that the 
mean field in the Hartree simulation is equal to the field 
in the classical simulation. Hartree simulations at small 
and large coupling will also be compared. In all cases 



^ Note that it is not possible to do this at the level of the action 
or Hamiltonian. 
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N=512, Lm=32, modes, E/Lm^=[0.01 :0.19], X/m^=1/1.25, u=0.0 

15 F ' ' ' ' ' ' ' 



E 




-15 I- , , ; 

2000 4000 6000 8000 10000 12000 14000 16000 
tm 

(a) Classical, Xjm?' = 1/1.25. For clarity the x-coordinate is 
shifted over L/2 through the periodic boundary conditions. 



N=512, Lm=32, 512 modes, E/Lm^=[0.01 :0.1], Vm^=1/1 .25, u=0.0 




(b) Hartree, A/m^ = 1/1.25. 




FIG. 3: Ener gy de nsity contour plot for kink-antikink annihilation from rest, for different couplings, classical and Hartree. All 
plots, except |3(c)| show total energy density. The vertical energy density range is chosen such as to emphasize the deviation 
from vacuum which in some of the plots will cut-off the peak at the center of the kink. 



the same physical volume will be used. Finally the effect 
of damping on the lifetime of the kink-antikink pair is 
discussed. The combined results are shown in Fig. 13 in 
the form of contour plots of the energy density for the 
4 different simulations. We have made plots of the en- 
ergy density instead of the (mean) field, as they show the 
position of the kink more clearly. 



Fig. 3(a) shows a classica l simu lation at the strong cou- 
pling Xjm^ = 1/1.25. Fig. |3(b)| sh ows t he same simula- 
tion, using Hartree dynamics. Fig. 3(c) shows the same 
Hartree simulation, with a damping T = 0.4to switched 
on during tm < 15. Finally Fig. |3(d)] shows the result for 
a Hartree simulation at a weaker coupling A/m^ = 1/12, 
without a damping term in the equations of motion. 



A comparison of classical and Hartree dynamics, 
3(b)| respectively, shows an enormous dif- 



Figs. 3(a 



and 



ference in annihilation time between them. The high 
instability of the Hartree kink antikink pair is caused by 
the fact that at this strong coupling, the average energy 
density e/m^ — 0.025 is higher than the potential barrier 
Ae/m^ — 0.019 between the three local minima (remem- 
ber that the Hartree approximation incorrectly predicts a 
first order phase transition instead of a crossover) , so by 
exciting particles and spreading the energy the barrier 
effectively disappears. The fact that we do not use an 
actual stationary solution to the equations of motion Q 
therefore causes the KK pair to evaporate very quickly. 
It is important to mention that, although it seems from 
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Fig. |3(b)| that the kink and antikink each spht into two, 
this is not actuaUy the case. By making an animation of 
the mean field, we found that it just "contracts" over the 
barrier, thereby annihilating both kink and antikink. In 
this process some energy is send off in opposite directions, 
but this cannot be seen as kink- antikink pairs. 

It is also interesting to note that forming a Bose Ein- 
stein distribution with the same average energy density 
would have an inverse temperature [3m — 3.5, just below 
the phase transition at which the three minima become 
degenerate. Combined with the low barrier this makes 
it profitable for the configuration to disintegrate, in or- 
der to lower the gradient energy, and settle down in the 
symmetric minimum. 



Comparing the Hartree results in Figs. 3(b) and 3(c) 



confirms the above interpretation. By damping the equa- 
tions of motion we establish two things. We make the 
configuration closer to a truly stationary one, thus mak- 
ing it more difficult to excite particles, and by lowering 
the average energy density the corresponding Bose Ein- 
stein temperature moves away from the phase transition 
point. 

After the damping is switched off at tm = 15, the kink- 
antikink pair survives until tm « 450. Switching off the 
damping at a later instant increases the survival time: 
using tm — 40 instead of tm — lb results in a survival 
time of tm « 500, i.e., the survival time increases by 50 
due to an increase in "switch-off time" of 25. When we 
do not switch off the damping, the pair survives until 
at least tm = 2000. We have not simulated this system 
longer, but we do not think it will ever decay. 

We also checked the influence of the precise value of the 
damping constant F, by making additional simulations at 
F = 0.2m and F = 0.8m. At the smaller value, F = 0.2m, 
the pair survival time decreases slightly, Aim = —20, at 
a larger damping, F = 0.8m, the survival time increases 
by only Aim — +2. Apparently, using a damping con- 
stant F = 0.2m until tm — 15 already removes most 
of the fluctuations around the stationary configuration. 
However, the remaining configuration is still much less 
stable than its classical counterpart. This is confirmed 
by a comparison of the total energy density with the de- 
picted mean field energy density: it turns out that the 
total energy density is much smoother than that of the 
mean field. Apparently the modes in part cancel out the 
kink-antikink inhomogeneities, thus facilitating the final 
annihilation of the pair. 

There is an important subtlety with the implementa- 
tion of the damping term. Since it is far from easy to 
add such a term to the equations of motion of the mode 
functions - even in the vacuum their unrenormalized ki- 
netic term is not vanishing ~ we only implement one in 
the mean field equation of motion. However, this makes 
the actual decrease of the energy density much slower, as 
it involves the transfer of energy from the modes to the 
mean field. We have seen in the simulations that there 
are basically two timescales, a rapid initial decrease, ap- 
proximately with the damping rate expected, after which 



most of the kinetic energy of the mean field has disap- 
peared, followed by a very slow decrease, draining the 
energy of the modes. When switching off the damping 
after the initial stage, as we have done in the simulation 
of Fig. 3(c) the kinetic energy of the mean field starts 
growing again, be it slowly, and a new energy equilib- 
rium is established. When the damping is not switched 
off, the system evolves exponentially slowly towards its 
final state. We will come back to this in the next section, 
in the study of the Hartree kink mass. 

,2 



By 



comparing Fig. 

,2 



3(b)1 at A/m^ = 1/1.25 with 



Fig. 3(d) at A/m — 1/12, we see that at the weaker cou- 
pling, the use of the classical kink form for the mean field 
and plane waves for the modes is much closer to a station- 
ary solution, even without damping. The configuration 
is oscillating around an approximately stable solution for 
an extensive period of time, before annihilation. By com- 
paring Fig. 3(d) with the damped strong coupling result 
Fig. 3(c) we see that the kink is broader at the stronger 



coupling, due to a change in the effective potential caused 
by the quantum modes (we also checked this broadening 
directly in a plot of the mean field) . Classically the width 
only depends on the mass m, which by construction is the 
same for both couplings. This change of width gives an- 
other interpretation for the fact that at larger coupling, 
the kink-antikink configuration is less stable: relative to 
their own size, the separation is smaller. 

As a final remark, it is important to note that the 
annihilation in Figs. 3(b) 3(d)| is not due to tunnelling 
but is an "over the barrier" annihilation. In Fig. 3(a) 
this is clear as it is a purely classical simulation, but 
one could suspect that the much shorter decay time in 
the quantum Hartree simulations is due to tunnelling. 
However, in all cases, including the classical, the decay 
is due to the attractive force between the kink and the 
antikink. The attractive potential falls off exponentially 
with the separation, and small deviations from the sta- 
tionary state only lead to decay after very long times. 
At strong coupling the deviations are larger and, as we 
have seen, in the Hartree simulations the effective rela- 
tive distance becomes smaller due to the broader kinks. 
In all runs in Figs.|3| except Fig. |3(b)1 the kink and an- 
tikink are moving towards each other before the decay, 
due to this attractive force, until they bounce and anni- 
hilate. After that the mean field e nds up in one of the 
In Fig 



3lb^ 



the barrier is so low 



two "broken minima" . 
that the field just moves over it and decays in that way. 
In this run, the field oscillates around the "symmetric 
minimum" after the decay, as already discussed before. 



C. Numerical results: Kink mass 

In this subsection we will discuss the kink mass, the 
energy of the static solitonic configuration. Classically it 
is given by Eq. H15|l . In the full theory this expression 
needs corrections. The first corrections were found by 
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Dashen et al.0|, see also Rajaraman j^,^ and more re- 
cently by Alonso Izquierdo et al.^^. A numerical study 
was done in Ref. 16]. The result is obtained using a 
semiclassical perturbative expansion around the classical 
kink solution. The kink mass is then given by the lowest 
energy level. This mass has to be renormalized and the 
net result is 

„,,,.2^ + „,,(^-i.)+e,(A), (19) 

A partial calculation of the order A correction can be 
found in Re f. [l7l . In this reference and independently 
also in Ref. 1^ a full calculation is done for the order A 
correction of the kink mass in the sine-Gordon model. 

More recently lattice Monte Carlo studies of the kink 
mass have been carried out [T^ |20| . In these papers 
two methods of calculating the kink mass were used: 
first using the fluctuation-fluctuation 2-point function 
{fj,{t) fi{0)) introduced in Refs. 0, which in imag- 
inary time should decay as exp(— Mkinki)- The second 
uses the difference in the ground state energy when peri- 
odic or antiperiodic boundary conditions are used. Both 
studies focus on a range of parameters, where Xa^ and 
—HqO? are of order unity, i.e., at large lattice distances. 

Very recently, Bergner and Bettencourt studied 
the kink mass in the Hartree approximation. Their equa- 
tions are derived from a 2PI effective action, but applying 
the Hartree approximation to it results in the same equa- 
tions of motion we are considering here. Using a numer- 
ical relaxation method they have obtained the Hartree 
kink mass at \/m? — 0.5 to be 0.4691 (translated in our 
units), which is about 70% of the classical value. We have 
reproduced their value using antiperiodic boundary con- 
ditions and a damping term ifT^ in the mean field equa- 
tion of motion. On a lattice with TV = 64 and Lm = 16, 
a/ao = 10 using a damping coefficient V /m = 0.4 we find 
at a time tm = 9 x 10^ - 10 x 10^ a value Mkink/'Ti = 0.47. 
Doubling the lattice size or decreasing the temporal lat- 
tice distance increases this value by about 1%, while in- 
creasing the simulation time with an order of magnitude 
decreases it by about 2%. Unlike Bergner and Betten- 
court we have not fixed the center of the kink, but the 
fact that the mass turns out equal shows that the zero 
mode can have at most a small effect, in agreement with 
Ref. [23. 

In order to make a detailed comparison with the full 
theory, i.e., the Monte Carlo data of Refs. JiL, J^l^i 
have done simulations for a range of parameters. The 



Note that there are two errors in 13: In Eq. (5.69) ^p^ + + 
4), in the numerator of the first term in the integral, should 
be (p^ + 1) Vp'^ "1-4 and \/p'^ + 2\ii the numerator of the second 
term in the integral should be ly^^ -|- 4. This second error comes 
from the transition of 2k^ — > . 
^ Note that the horizontal axis, A, in Fig. 4 of Ref. |20| should be 
multiplied by 6, as can be derived from the curve of the semi- 
classical kink mass. 



N=64, ao/a=0.1, r/m=0.4, tm=90x10^... 100x10^ 




(a) a'^^? = -1.0. 



N=64, ao/a=0.1, r/m=0.4, tm=90x1o'... 100x10^ 




(b) o?^? = -2.2. 



N=64, ao/a=0.1, r/m=0.4, tm=90x10^... 100x10^ 

4 I ■ 

Hartree, m.f. a ji|.gn=-4.0 — 

Classical — 

3 - Hartree, — 

Semi-classical 




3 3.q 4 4.5 5 5.5 



(c) a2^2 ^ _4 Q 



FIG. 4: Comparison of classical, semiclassical (one loop) and 
Hartree results. The Hartree results are plotted for constant 
renormalized /.tren and unrenormalized ^Mq. Also see the text. 
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results are shown in Fig. 2| We show results for constant 
renormalized fj,^^^-^ as well as constant unrenormalized /Iq. 
For the constant /i^g„ we also show separately the mean 
field contribution, defined as in Eq. (16), Ref. 0. The 
data in Refs. [l^ l20l | is shown as a function of the un- 
renormalized /ip, in order to study the renormalization 
effects. However, to make a meaningful comparison with 
the classical and semiclassical 1-loop results, it seems 
more proper to compare at constant yLt^en- The arrows 
in the plots show the zero temperature phase transition 
point, at which the three minima become degenerate, cal- 
culated using the static effective potential, cf. Ref. 
The dots indicate the point at which the "broken" min- 
ima vanish altogether. In order to make the compari- 
son with the Monte Carlo data easier, lattice units have 
been used. Most simulations where done with N = 64, 
ao/a — 0.1 with a damping term F/m = 0.4 while the en- 
ergy was measured at tm « 100 x 10"^. We have checked 
for lattice artefacts (the influence of N and ag/a), their 
combined error is not larger than about 1%. Also the 
effects of the precise value of the damping term and time 
at which the mass is measured are small and of the same 
order of magnitude as that of lattice artefacts. 

We see that for all renormalized masses, the Hartree 
results are between the classical and semiclassical values. 
Using the unrenormalized mass as the parameter to com- 
pare with, the Hartree results are much lower than the 
semiclassical values. Comparing these with Figs. 4-6 of 
Ref. we find that the Monte Carlo results are above 
the Hartree results. Note that the semiclassical line in 
Fig. 6 of Ref. ^| is approximately 0.3 too high. For the 
two higher values of |/ironl j where the data taken is further 
away from the first order phase transition, the mean field 
contribution is almost equal to the classical value. The 
difference between the Hartree and classical kink mass 
only comes from the modes. For the lower a^|/i^Q„| = 1.0 
the "mean field mass" behaves very interestingly. We see 
that, coming from below, it goes up again near the phase 
transition. This is caused by the emergence of the sym- 
metric minimum. More data would be needed to study 
the precise behaviour at the phase transition, but from 
its definition, Eq. (16a) in Ref. it seems it only goes 
up a final amount, due to this symmetric minimum and 
the change in the effective potential. It is also not clear 
to us at this stage whether it is only an artefact of the 
way we split the total energy in a mean field and mode 
contribution or is an actual physical phenomenon. 

Finally, we have measured the mass at a coupling 
X/rn^ = 1/1.25, used for example in Figs. 3(b) and |3(c)| 
using Lm = 16, = 64 at a time tm = 2 x 10 '^. In prin- 
ciple we could find the value also in Fig. |4(a)| but those 



data where obtained a large lattice distance and a sepa- 
rate test is useful. The result is Mkink/"T. — 0.22, or 54% 
of the classical value .42. This value for the kink mass 
is consistent with Fig. |4(a)| showing that we are indeed 
in the scaling region, as claimed by Ref. [23| (the same 
is also found for the weaker coupling Xjrr? = 0.5, con- 
sidered in Ref. [l^). The semiclassical value (|19|) from 



Ref. [Ijl gives a value 0.08, about 1/5 of the classical 
value. Note that at X/m^ ~ 1 the semiclassical approxi- 
mation breaks down, showing up by a change of sign in 
the kink mass, the order A'' correction becomes equal to 
the leading term. In the Hartree approximation we 
can in principle continue till A/m^ w 1.2, where the zero 
temperature phase transition occurs. 



IV. MOVING HARTREE KINKS 
A. Initial condition 

In order to describe kink-antikink collisions, it is nec- 
essary to have a description of a moving Hartree kink as 
well. In the full quantum theory this poses a difficult 
problem, since we do not know how to Lorentz trans- 
form the "quantum cloud" surrounding the kink. In the 
Hartree approximation we are saved by the fact that the 
field is completely expressed in terms of ordinary func- 
tions of X and i, and we can find a moving kink by per- 
forming a simple Lorentz transformation: 

fix, t) f'(x, t) = f{^[x - ut],^[t ~ ux]), (20) 

both on the mode functions and the mean field. 

However, there are still a number of difficulties with 
finding a proper initial condition for a kink-antikink pair, 
moving toward each other. First of all, from Eq. H20|l it 
follows that we need the mode functions on a backwards 
space-time line, while we do not know them in analytic 
form, only from numerical simulations. For the classical 
kink solution this problem does not exist, since it is not 
only stationary, but also static and known analytically. 
Second, we have to boost the kink and antikink, with 
their respective mode-functions, in opposite directions. 
This means that they have to be combined in a nontriv- 
ial way in the middle, they "shift into each other" . It also 
means that we have to double the density of mode func- 
tions in fc-space, since the physical space doubles (the 
best one can do is determine the mode functions for a 
single kink). We have to decide how to do this in a con- 
sistent way. Finally there is the problem that the boosted 
coordinates ^{x — ut) and 7(t — ux) will in general not 
fall onto space-time lattice points in the unboosted frame 
or vice versa. 

In order to circumvent these problems, we will just 
use the vacuum form as initial conditions for the mode 
functions. Since the vacuum is invariant under Lorentz 
transformations, we can use the same mode functions in a 
boosted frame. As mentioned before, the mode functions 
will be close to the vacuum form, if the vacuum expecta- 
tion value V is large. Furthermore by boosting the kink 
solutions we reduce the width of the kinks. However, in 
the center of mass frame, we keep the lattice size the 
same thus effectively increasing the relative distance be- 
tween them. This leads to much more stable solutions, 
which is even further enhanced by the time dilatation. 



9 



Before we proceed we will outline a set of solutions 
for the aforementioned problems, which will yield a very 
accurate, boosted, Hartree kink-antikink solution and 
will therefore also be usable at strong coupling and low 
speeds. In the simulations we will in general, for simplic- 
ity, only use the vacuum form mode functions. 

The way to solve the stability problem of a (single) 
kink is to use damped equations of motion, as we already 
used in the preceding subsections, especially in combina- 
tion with antiperiodic boundary conditions on the mean 
field to make the kink absolutely stable. The damping 
will give us stationary, but time-dependent solutions for 
the mode functions in the background of a kink, which 
subsequently can be Lorentz boosted. 

The problem of the boosted coordinates not resulting 
from discrete lattice points in the unboosted coordinate 
can be adequately solved by linear interpolation, pro- 
vided the lattice distance is not too large. The related 
matching problem is not a real problem: we have to 
match the two functions at two consecutive time steps 
and when they approach each other with speed u, they 
shift into each other over a distance uao <C a, i.e., the 
mismatch is less than a lattice distance and the linear 
interpolation automatically solves the problem. 

The remaining problem is to combine the two sets of 
mode functions, of the separate kink and antikink config- 
urations, into one set which is twice as large. One possi- 
ble solution is to determine both kink and antikink in a 
volume that already has the size of the combined config- 
uration and then use an averaging procedure to combine 
the modes. This has the disadvantage that oscillatory 
functions with different phases are added. Another way, 
which we suggest, is to determine the two configurations 
in their original volume and combine them afterwards 
into one set in the double volume. When using periodic 
mode functions, it follows from the antiperiodicity of the 
mean field that at all times 



fk{x,t) = f-k{L - x,t), 



(21) 



provided this relation holds aA, t — and t — ao, a condi- 
tion which is true for the plane wave initial conditions we 
use. Note that at later times the k label is no longer equal 
to the Fourier label, for which the relation holds trivially. 
We can therefore combine each mode function fk with 
either itself or with f-k, since both combinations yield 
continuous functions. The derivative will in general be 
discontinuous by changing sign, possibly causing trouble 
in the equation of motion. However, the energy density 
will be finite and even continuous, as it only depends on 
the square of the first derivative, which is continuous. 

We end this subsection summarizing the above steps 
one could follow to obtain a better initial state for a 
colliding kink-antikink pair. First, using antiperiodic 
boundary conditions for the mean field and a damping 
term —Tdt4>, o^s evolves the Hartree equations, start- 
ing from a classical kink with free field mode functions. 
One thus obtains a stationary Hartree solution for a sin- 
gle kink, but has to keep track of the modes and mean 



field on a backward space-time line. Then one constructs 
the mean field from the kink and its mirror image. The 
modes are constructed by combining fk{x,t) with itself 
and with f^k{x, t) — fk{L — x, t), thereby doubling their 
number. 



fk,i(x,t) 

fk,2{x,t) 



fkix,t), 0<x<L, 
Jk{x-L,t), L<x<2L, 

' fk{x,t), 0<x<L, 
f-k(x-L,t), L<x<2L. 



(22a) 
(22b) 



The latter, symmetric, combination can lead to functions 
in which the second derivative diverges as 1 /a: the energy 
density, depending only on the square of the first deriva- 
tive, is finite and continuous in the continuum limit, how- 
ever it might cause problems in the equations of motion. 



B. Numerical results: Thermalization from 
colliding kink-antikinks 

We start the discussion of our numerical results with 
the thermalization properties of a colliding kink-antikink 
pair, using the defining relation for the particle num- 
ber Uk- We used a small coupling, X/m'^ = 1/12, in 
a volume Lm = 32 at a speed for which 7 = 2. The 
classical kink mass at this coupling is Mkink/w = 4. 
Since we initialize the mean field as a classical kink and 
the modes in the vacuum, the initial energy density is 
E/Lm^ = 2 • (4 + 4)/32 = 0.5. At this energy density 
and coupling we found, in Ref. both in dynamical 
simulations and an effective potential calculation, a tem- 
perature T/m « 1.0. 

In order to enable us to correctly interpret the result, 
we also measured the "particle number" at tm = 0. Of 
course at this time the system is so far from equilibrium 
that we cannot interpret the result as quasiparticles, but 
it does give an idea of the energy distribution as a func- 
tion of k. The result is plotted in Fig. |S1 The similar- 
ity with a Bose-Einstcin distribution is remarkable and 
we have to be careful in interpreting the results at later 
times. One of the indications, that it is not a truly ther- 
mal distribution, is its very low "temperature": at an 
energy density E/Lrri^ = 0.5 we expect to find a temper- 
ature around T/m = 1.0, while Fig.|31seems to indicate a 
value which is almost a factor of 3 smaller. Furthermore, 
the particle number mainly comes from the even modes, 
as a result of the symmetric initial condition. 

Fig.0shows the actual particle number, at later times, 
when the kink-antikink have already annihilated. Then, 
we do find a temperature around 1, conform the results of 
Ref. As a comparison we also plotted, in the same fig- 
ure, the result obtained in that paper using a flat ensem- 
ble simulation. It clearly shows that the temperatures 
are equal: T is uniquely determined by the coupling and 
energy, not by the initial condition. The initial condition 
still has an influence in determining which modes have 
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T N=256, Lm=32, Vm^=1/12, E/Lm^=0.5 




FIG. 5; Initial form (tm — 0) of the particle distribution for 
a Hartree kink-antikink collision. The initial form already 
seems close to a Bose-Einstein distribution. 



T N=256, Lm=32, Vm^=1/12, tm=9000.. 10000, E/Lm^=0.5 




0.5 1 1.5 2 2.5 3 3.5 % ^ 



FIG. 6: Late time form of the particle distribution for a 
Hartree kink-antikink collision, tm ~ 10000, after annihi- 
lation. 



thermalized. We see that starting from an annihilating 
KK actually allows the system to thermalize faster, as 
it has a more favourable energy distribution, cf. Fig 
However, the KK initial condition has a delta function 
initial density matrix, which is not very suitable for our 
Hartree ensemble approximation: we do not reduce the 
statistical errors by averaging over multiple initial condi- 
tions. 

Having showed that this initial condition also leads 
to an approximately thermal Bose-Einstein distribution 
after annihilation, we will now leave the thermalization 
topic and continue with a discussion of the actual colli- 
sions. 



C. Numerical results: Critical speed 

We start this section with a study of classical kink- 
antikink collisions. An extensive study of this can be 
found in Ref. 24]. 

Figs. [7| show energy contour plots of Lorentz boosted 
classical kink-antikink configurations. The results in 
these plots demonstrate how, at very low incident speed, 
they just annihilate, while at a certain range of higher 
speeds, they bounce a few times and then either annihi- 
late or escape again to infinity, at a slightly lower speed. 
At high enough speeds, when u is larger than a certain 
Uc, which has a value between 0.25 and 0.30, the pair im- 
mediately escapes after the collision, again with a lower 
speed. Part of the kinetic energy is transferred to an in- 
ternal vibrational mode of the kink, which can be seen 
as a wiggling in the contour plots. 

These results are consistent with Ref. |2^ in which 
a critical speed Uc = 0.2598 is found, above which col- 
liding kinks always bounce back, while below a speed 
0.193, they always annihilate immediately. Between 
these speeds Ref. found reflection and annihilation 
bands, caused by a resonance between the center of mass 
motion and the vibrational mode. This was confirmed 
in Ref. 25], also showing that, when zooming in on the 
bands, their structure behaves fractally. This behaviour 
could even be reproduced using a collective coordinates 
approach, reducing the infinite-dimensional system to a 
two dimensional Hamiltonian one. 

Except for Fig.| 7(d)]all shown results are in accordance 
with Table I in Ref. |24l]. Since we used different dis- 
cretizations, etc., and given the very narrow width of the 
stability bands already found by Ref. (2^] and the frac- 
tal behaviour found in Ref. J^^, we can certainly expect 
small differences in the precise position of these bands, 
explaining the discrepancy. The authors of Ref. 
found that the final speed, for an initial speed above Uc, 
satisfies the following relation 

uj (X uj — u^, (23) 

which is consistent with Fig. |7(e)| from which we derive 
a final speed Uf w 0.13 — 0.14. For Ui = 0.30 relation 
(ESI predicts Uf = 0.15. 

It is important to note that although breatherlike 
states seem to emerge after a collision, truly stable 
breathers do not exist in the theory |20]. However, 
very long-lived and almost stationary configurations do 
exist [23 as we have shown here again. The radiation 
causes the energy to decay as l/log(i), as found analyti- 
cally by Ref. 26] and confirmed numerically by Ref. ^3] ■ 

We would like to compare these classical results with 
the Hartree approximation to see if, in the quantum the- 
ory, a critical speed still exists. In the classical theory, 
one can always use units such that /i = A = 1, and one 
expects a unique critical speed, while in the quantum 
theory, including the Hartree approximation, this is not 
the case. However, we will not go much further into the 
question of the coupling dependence of the critical speed. 
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N=512, Lm=64, modes, E/Lm =[0.085:1 .7], X/m =1/12, u=0.19 



N=512, Lm=64, modes, E/Lm =[0.085:1.7], X/m =1/12, u=0.20 
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(a) Ui = 0.19 (7 ss 1.019). 



(b) Ui = 0.20 (7 fs 1.021). 



N=512, Lm=64, modes, E/Lm =[0.085:1 .7], X/m =1/12, u=0.21 



N=512, Lm=64, modes, E/Lm =[0.09:1 .8], X/m^=1/12, u=0.25 
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(c) Ui = 0.21 (7 ^ 1.023). 



(d) Ui = 0.25 (7 fs 1.033). 



N=512, Lm=64, modes, E/Lm =[0.1 :2], X/m^=1/12, u=0.30 




(e) Ui = 0.30 (7 « 1.048). 



FIG. 7: Colliding classical kink-antikink, A/m^ = 1/12, above and below the critical speed. 
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N=512, Lm=64, 512 modes, E/Lm =[0.15:3], X/m =1/12, u=0.70 



N=512, Lm=64, 512 modes, E/Lm =[0.175:3.5], X/m =1/12, u=0.75 




(a) Ui = 0.700 (7 » 1.400). 



(b) Ui = 0.750 (7 « 1.512). 



N=512, Lm=64, 512 modes, E/Lm =[0.15:3.5], XVm =1/12, u=0.76 



N=512, Lm=64, 512 modes, E/Lm =[0.15:3.5], X/m =1/12, u=0.765 
30 




(c) Ui = 0.760 (7 Ri 1.539). 



(d) Ui = 0.765 (7 ^ 1.553). 



N=512, Lm=64, 512 modes, E/Lm =[0.15:3.5], X/m =1/12, u=0.77 



N=512, Lm=64, 512 modes, E/Lm =[0.2:4], X/m =1/12, u=0.80 




(e) Ui = 0.770 (7 » 1.567). 



(f) Ui = 0.800 (7 » 1.667). 



FIG. 8: Colliding Hartree kinli-antikinli, X/m^ = 1/12, above and below the critical speed. 
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N=512, Lm=64, Wm =1/12, Uc=0.7607 
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FIG. 9: u^f as a function of uf in the Hartree approximation, 
showing roughly the same behaviour as in classical theory, 
but with a itc = 0.7607. 



In Figs. IHl the Hartree results for six different initial 
speeds are shown. We see that a critical speed does in- 
deed exist, with a value somewhere between 0.760 and 
0.765, indicating that the quantum kink pair is less sta- 
ble than its classical counterpart, at least at this coupling. 
One might think that this instability is caused by the fact 
that we only have an approximate quantum kink-antikink 
as initial condition, but as can be seen from Figs. |H1 the 
pair before the collision looks very stable, there is hardly 
any wiggling, while after the collision the wave packet 
is dispersing and oscillating and its speed has decreased 
considerably. The reason for the higher instability is the 
radiative channel. This probably also causes the fewer 
number of breather states found: only at an initial speed 
u — 0.760, very close to the critical speed, can we recog- 
nise an approximate breather state. Bands of stability 
have not been found and in most simulations in which 
the pair does not bounce back immediately, the radia- 
tion causes them to annihilate. 

We also have investigated the functional behaviour of 
as a function of u^, to see if it has the same form (|23ll . 
valid in the classical approximation. The result of several 
runs at different initial speeds is shown in Fig. O We see 
that close to the critical speed, Uf does indeed behave in 
the same way as in the classical theory. 



5.34 (i 



(24) 



The resulting critical speed Uc can be determined very 
precisely in this way and we obtain Uc = 0.760. Further 
away from the critical speed the behaviour is not linear in 
uf, in contrast to the classical result, but the corrections 
are relatively small. 

We have done a few simulations using a stronger cou- 
pling A/m^ = 1/6 and they indicate the general be- 
haviour is only slightly different. For example, the critical 
speed is not very different, it has a value between 0.790 
and 0.795, from a linear fit of uj versus uf we obtain 



Uc = 0.793 instead of the value 0.7607 at small coupling. 
Furthermore at a Ui = 0.80 we find Uf ~ 0.280, while at 
A/to^ = 1/12 we found Uf = 0.538. Since Ui = 0.80 is 
much closer to this value of Uc, it explains why the value 
oi Uf a,t Ui — 0.80 is a factor of 2 smaller. The prefactor 
in Eq. H24|l is higher at the stronger coupling, around 7 to 
7.5 instead of the 5.34, but more data would be needed 
in order to obtain an accurate answer. 

Apart from the factor 3 larger value of the critical 
speed, their are other important differences between the 
classical and Hartree results. Figs. |H| and |7| For exam- 
ple, irrespectively if the pair annihilates or not, a lot of 
energy is radiated away after the collision in the form 
of quasiparticles, moving with approximately the speed 
of light. One might think that the radiation is mainly 
described by the mode function contribution to the en- 
ergy density, while the energy density of the surviving 
kink pair comes from the mean field. However, by check- 
ing the different contributions separately we found this 
is not the case, the kinks are mainly described by the 
mean field, but also have a contribution from the modes 
and the radiation is described by both together. Only 
the total field is a physical quantity describing both the 
kink-antikink and the radiated particles. Just as only 
the total two point functions describe the quasiparticles 
which approximately thermalize. 



D. Numerical results: Scaling 

As we can describe kink-antikink collisions in a quan- 
tum theory, it is interesting to study the possibility of 
using them as a toy model for heavy ion collisions, such 
as has been conducted at the SPS, are currently being 
conducted at RHIC and will later be at the LHC. We 
therefore look closer at the region just after the actual 
impact, in collisions with high 7. In Figs.^jwe show the 
results of four different simulations, at different 7 factors 
and couplings, including one using classical dynamics. 

The plots show a remarkable similarity. In order to 
investigate this similarity more quantitatively we made 
time slices shortly after the impact, shown in Figs. 1111 
at a time 11.5, the impact itself was at a time tm « 8 
(the determination of the exact time of impact has some 
complications to which we will return later). 

In Fig. 



11(a) we see that a different impact velocity 



only influences the resulting energy density in the emerg- 
ing kink-antikink, the central region is the same for both 
speeds. There, the difference in initial speed only shows 
up in the oscillation wavelength, which scales due to the 
7 factor. 

In Fig. ^] we have plotted a similar comparison at 
smaller coupling, X/m? = 1/12, as a function of space- 
time rapidity rj — arctanh(z /t) at constant "proper time" 
— lOm^^ {t — at the collision). Even on 



the linear vertical scale, the plateau heights arc almost 
identical. It is also striking to see how flat the plateau 
is, an indication for hydrodynamic scaling. Note that 
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N=1 024, Lm=32, 1 024 modes, E/Lm^=[0.01 :0.52], X/m^=^n .25, Y=20 
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N=1024, Lm=32, modes, E/Lm^=[0.01 :0.52], X/m^=^n .25, 7=20 




(a) Hartree, A/m^ = 1/1.25, 7 = 20 (« = 0.9987). 



(b) Classical, A/m^ = 1/1.25, 7 = 20 (u = 0.9987). 



N=51 2, Lm=32, 51 2 modes, E/Lm^=[0.01 :0.52], X/m^=1 /1 .25, y=1 



N=512, Lm=32, 512 modes, E/Lm^=[0.1 :5], X/m^=1/12, Y=10 
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(c) Hartree, A/m^ = 1/1.25, 7 = 10 (m = 0.9950). 



(d) Hartree, A/m^ = 1/12, 7 = 10 (« = 0.9950). 



FIG. 10: Colliding kink-antikink, period around the collision only. 



xm « 10 sinh(l) « 12 for the edges of the plateau, a 
large distance compared to the size of the kink, which is 
much smaller than 1. 



In Fig. 11(b) we show the difference between the two 
couplings. We have rescaled the energy density with 
its average over the whole system, in order to facili- 
tate the comparison. Again, the difference is small. At 
the stronger coupling the energy is slightly more concen- 
trated in the kink-antikink. 



Finally in Fig. 11(c) we compare the difference between 
Hartree and classical dynamics. In both simulations the 
kink-antikink pair region is very similar. The central 
region shows some differences, although the total Hartree 
energy density in this region is very close to the classical 
energy density. The modes are essential for this, the 
mean field is more concentrated in the kink-antikink. 

To summarize, we have shown that the collisions show 
scaling behaviour for an extensive range of energies and 



couplings: they result in approximately the same central 
region, which is very flat in space-time rapidity. Fur- 
thermore at high incident speeds the difference between 
Hartree and classical is relatively small: classical dynam- 
ics can be used in the study of these collisions, giving 
foundation for its use in the study of heavy ion collisions. 

It is interesting to compare the results obtained here in 
kink-antikink collisions with what is known about heavy 
ion collisions. It is already long accepted that the 
height of the central plateau only weakly depends on the 
incident speed. Experiments at RHIC, at center of mass 
energies of 130 and 200 GeV per nucleus (i.e., 7 — 130 
and 7 = 200, respectively), indicate an increase in par- 
ticle density of 14 ± 5% [23 • At the smaller coupling of 
A/m^ = 1/12 we find, at a time tm = 3 after the colli- 
sion, an increase of only about 3% when increasing from 
7 = 10 to 7 = 20. However, given the very different 
physical systems, a better match would be a mere coinci- 
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(a) Hartree kink-antikink at 7 = 20 and 10 
{X/m? = 1/1.25). 
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FIG. 12: Total energy density in "Bjorken variables" r and 77 
for two different initial speeds. Note the linear scale. 
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dence. Longer after the collision, mr ~ 10, ... , 100, the 
difference in energy density between the two simulations 
becomes somewhat larger, as we will see later. 

Furthermore, our findings, that only a relatively small 
fraction of the energy remain in the central region and 
the concept of a central plateau in space-time rapidity 
are, at least qualitatively, consistent with real-life heavy 
ion collisions. A quantitative comparison would not be 
very meaningful. 

A first study of hydrodynamic scaling in a 0"* theory 
can be found in Ref. [20|- The authors calculate the 
energy momentum tensor, which in our metric equals 



(b) Hartree kink-antikink at \/vn? = 1/12 and 1/1.25 
(7 = 10). The energy density is normalised with the 
average energy density. 



I N=1024, Lm=32, am=1/32, A;m^=1/1.25, tm=11.5, y=20 



10 




Hartree, mean field 
Hartree, total 
Classical 




(c) Hartree mean field, Hartree total and classical 
(A/m^ = 1/1.25, 7 = 20). 



(25) 



in two systems, a colliding kink-antikink^ and a decaying 
Gaussian wave packet. For a perfect fluid the energy 
momentum tensor can be expressed in the energy density 
and pressure. The assumption of a perfect fluid is valid 
when collisions can be neglected, as in the (homogeneous) 
Hartree approximation. For example, the trace of T^i, 
gives (in 1 -f 1 dimensions) 



n = -e + p=-2V, 



1 



(26a) 
(26b) 

(26c) 



Note that this is a slightly different definition from the 
one given in Ref. "so], Eq. (13), as we would like both Too 
and Til to vanish in the vacuum d> = v. Also note that 



.11: Comparison of energy density profiles at a time-slice 
— 11.5, (collision at tm — 8) for different parameters. 



Note that in Ref. l3fl a product of a kink and antikink is taken, 
while we use a sum of the two, Eq. I17i . Of course both are 
approximate kink-antikink solutions and the difference should 
be small. 
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(a) Log-log plot of energy density at a; = 0. 



(b) Pressure p = e — 2V and energy density. Note that e 
falls right on top of p. 



FIG. 13: Energy density and pressure in the central region, after a classical kink-antikink collision, 
consistent with p = constant and Co = 0. 



At later times, both are 
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FIG. 14: As in Fig. |13(a)| but now for a Hartree kink-antikink 
collision. 



we just write the classical expressions, the extension to 
Hartree is straightforward. From Eqs. (|26|l we find the 
pressure 



(27) 



Using the scaling behaviour one can then, for example, 
express the speed of sound in the energy density in the 
center 



e cx T 



-(1+-;^) _ at X = 0, 



(28) 



where r = \^t^ — x'^ is the "proper time" and cq is the 
speed of sound. For a classical kink-antikink collision at 
strong coupling A/m^ = 1/1.25 we plotted the central 



energy density as a function of time (where the time is 
set to zero at the end of the collision), on a log- log plot 
in Fig. 13(a) and, on a linear scale, together with the cor- 
responding pressure, in Fig. 13(b) Both from Eqs. 



and Eqs. (|28l) we find that shortly after the collision, the 
pressure in the central region becomes zero, leading to a 
vanishing speed of sound. This indicates that after this 
initial stage, the interaction becomes negligible and we 
just see free expansion, in agreement with the flat cen- 
tral plateau. Initially, however, the system is interacting 
and there is a non- vanishing speed of sound of about 0.8. 
The actual value depends quite strongly on the time ori- 
gin. We have derived the value to use, from contour plots 
such as Figs. They show that the kink-antikink stay 
together for a short period of time Aim w 0.5. 

The result obtained from Hartree dynamics is shown 
in Fig. El In this case we can only find a speed of sound 
in the first stage, resulting in a value very similar to the 
classical result. The second stage exhibits strange oscil- 
lations but no power behaviour. This is caused by the 
emergence of a "symmetric minimum" at cj) ~ 0, already 
discussed in Sec. IIII Bl In Figs. E| we show the results 
for collisions at the smaller coupling A/m^ — 1/12. Here 
we see that both classical and Hartree results give an ap- 
proximately vanishing speed of sound and thus pressure. 
The initial speed of sound is for both smaller, about a 
factor of 2 compared to the strong coupling, which is 
(qualitatively) consistent with a smaller coupling. The 
Hartree result. Fig. |15(bj] shows small deviations at the 
latest times, indicating some interaction. At a higher in- 
cident speed, 7 — 20, the prefactor in the 1/t decay rises 
to 4.4, both for the classical and Hartree simulations. It 
is interesting to see that this increase is actually very 
close to the RHIC results as discussed above. 



17 



am=0.0625, Wm'^=1/12, y=10, classical 



am=0.0625, Wm =1/12, y=10, Hartree (all modes) 
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(a) Classical. 



(b) Hartree. 



FIG. 15: As in Fig. |13(a)| andlTllbut now for A/m^ = 1/12. 



It is difficult to compare our results quantitatively with 
those of Bettencourt et al. 59| : since their result for the 
speed of sound was obtained for a disintegrating Gaus- 
sian wave packet in the "symmetric phase." The effective 
coupling in the "symmetric phase" is much smaller than 
in the "broken phase" j^] and we have studied only the 
central region in the wake of a kink-antikink collision, as 
this is more similar to a heavy ion collision. The free 
expanding region is not present in the results of Ref. [s^l 
and the power behaviour is less pronounced. However, 
the speed of sound in the initial region is very compara- 
ble. 



boundary conditions. Furthermore, in the quantum the- 
ory, one might think the kinks will be fully described by 
the mean field, while the mode functions just describe 
fluctuations around them. However, as we already no- 
ticed before, this split-up cannot be taken so rigorously: 
the kink number will also be partially described by the 
modes and a definition based on the ensemble descrip- 
tion of our system is needed, as only ensemble averaged 
quantities are physically meaningful. 

In the classical theory, a useful quantity to look at is 
the following. 



; class 



(29) 



V. THERMAL KINK NUCLEATION 

In this section we briefly discuss the thermal creation 
and annihilation properties of the kink-antikink pairs. 
There is a long history of papers on the subject of ther- 
mal nucleation. See, for example, Ref. 31] and refer- 
ences therein for an analytical study, and for example 
Refs. [s^, I33L Is^ I35I Is^ for numerical studies. Very re- 
cently Ref. [33 has made a study in a heavy ion collision 
context, calculating classically the multiplicities of KK 
in an expanding Bjorken frame, starting from a Boltz- 
mann distribution. 

All these (numerical) studies so far have only consid- 
ered the classical nucleation rate. The Hartree ensemble 
approximation method allows us to study the creation of 
kink-antikink pairs starting from a thermal Bosc-Einstein 
distribution and to compare the classical and Hartree 
approximation. One of the problems we are thus faced 
with is a proper definition of kink number. Since only 
pairs, with no net kink number can be created, we need 
an effective kink number definition. The actual wind- 
ing number will always be or 1, depending only on the 



where the overline denotes a spatial average over the full 
lattice size. When a kink-antikink pair is present, this 
observable will become of the order u^. If the coupling is 
not too small, it will therefore give a reasonable indica- 
tion for their presence. Note that it cannot distinguish 
between one or multiple pairs and that it becomes smaller 
if the pair is not at maximum separation. The advantage 
is that we can easily extend its definition to the quantum 
theory 



(30) 



and that it can be easily found from Sk, Eq. It is 
however UV divergent and should be renormalized by 
subtracting the vacuum contribution: 



1 



(31) 



We now have a quantum kink indicator, showing if a 
pair is present and which can be easily compared to 
the classical theory. This is done in two simulations at 
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|3m=0.563, Vm^=1/12, N=256, Lm=64 
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FIG. 16: Kink indicator parameter Q from Eq. 1311 for classi- 
cal dynamics (upper plot) and Hartree dynamics (lower plot). 
Both values are normalised with the value 0.875, the value of 
Q for a symmetrically placed kink-antikink pair in a volume 
Lm = 64. Of course this value heavily depends on the exact 
position of the kinks. 



\/m? = 1/12 in a volume Lm = 64, at an inverse tem- 
perature f3m = 0.563, just below the thermal phase tran- 
sition at /3to — 0.562. At such a high temperature we 
expect the highest creation rate. In a volume Lm = 64, 
it would follow from Refs. 0, [s^, that we should find 
about a half to two kinks per total volume, depending 
on the counting algorithm. The result for Q is shown in 
Fig. El We can clearly see it is different from zero, show- 
ing the presence of kink-antikink pairs. Furthermore, 
these pairs emerge through the dynamics, initially the 
number is lower. Finally, although the quantum number 
is smaller than the classical, it does not go away, while, 
if we look at the mean fields in different realisations sep- 
arately, we find the kink-antikinks do seem to disappear 
in the Hartree approximation. From the Q indicator we 
conclude that they first are contained in the mean field, 
while later their description is taken over by the modes, 
again showing that only the total two-point functions, 
from the total field, describe the physical quantum field. 
As an example of their emergence from the dynamics, we 
plotted, in Fig. El the mean field of one of the realisa- 
tions at 4 different times. Initially the field fluctuates 
around one of its minima, with a reasonable amount of 
energy, it then forms a KK pair, which subsequently an- 
nihilates, while transferring its energy to the quantum 
modes, as can be seen from the smaller amplitude of the 
fluctuations. The total two-point function still describes 
KK pairs, but this cannot be seen from the separate 
realisations. 



VI. CONCLUSION 

In this paper we looked at the topological defects which 
can be present in a scalar Xip^ theory: kinks and an- 
tikinks. After rederiving the classical solutions, we used 
both classical and Hartree dynamics in studying their an- 
nihilation, both when initially at rest and when boosted. 
We found that classical kink pairs are much more sta- 
ble than their Hartree counterparts. Hartree kink pairs 
at smaller coupling are therefore also more stable than 
at larger couplings. Although an exact Hartree solution 
to the equations of motion is not known, by damping 
the mean field equation, we are able to find an approx- 
imate solution. This damping can, especially at larger 
couplings, prolong the survival time enormously, but is 
still orders of magnitude shorter than for a classical kink- 
antikink pair, due to radiation in the form of quantum 
particles, the quasiparticles discussed in Refs. We 
suggested an algorithm to find an accurate, numerical 
and more stable solution for the Hartree equations of 
motion, but have not yet implemented it. 

Using a damped mean field equation of motion we were 
also able to determine approximately the quantum kink 
mass. The Hartree approximation seems to give better 
results than the one-loop semiclassical approximation, es- 
pecially at larger couplings, where the semiclassical ap- 
proximation becomes less reliable. The mean field contri- 
bution shows a funny increase as a function of coupling, 
when approaching the phase transition, although the to- 
tal kink mass behaves regular. It is not yet clear how 
to interpret this. Damping the equations is not very ef- 
ficient, as it can only be implemented in the mean field 
equation, leading to very long relaxation times. We found 
a rapid short relaxation until a time tm w 10 — 20, in 
which most of the mean field energy is removed, followed 
by an exponentially slow relaxation in which the energy 
of the modes is drained indirectly via the mean field. 
Comparing the Hartree results with Monte Carlo results 
from Refs. 19, .2Qj seems to indicate that the Hartree 
results, although larger than the one-loop semiclassical 
results, are still too small. 

We have shown that a colliding kink-antikink pair, af- 
ter annihilation, also leads to an approximate thermal 
Bose-Einstein spectrum, just as a flat initial ensemble. 
The initial energy distribution has a form which ap- 
proaches the thermal distribution more easily than the 
fiat ensemble, as used in Ref. 0. This makes it possi- 
ble to recognise the Bose-Einstein distribution already in 
an early stage. However the statistical errors are larger, 
since we cannot average over multiple initial conditions. 
It does show that the result of Ref. Q is rigorous, it does 
not depend on the specifics of the initial ensemble. 

In the study of kink collisions, we found once more that 
the classical kinks are much more stable: for the classical 
dynamics, we reproduced the results of Ref. on the 
critical speed and the existence of approximate breather 
modes. In the Hartree approximation the critical speed is 
considerably higher and we were not able to find stability 
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FIG. 17: Thermal kink nucleation and subsequent annihilation in the Hartree approximation. The dotted lines are at the zero 
temperature "broken" and "symmetric minima" ±v{T — 0) and 0. 



bands, i.e., approximate breather solutions. In order to 
make this result rigorous, many more simulations have 
to be done, at higher numerical precision, since the ra- 
diation makes it difficult to see clearly if a short-lived 
bound state, or breather, has emerged. The precise de- 
pendence on the coupling is also still an open question. 
In the classical theory the dimensionless critical speed is 
independent of A as the action can be rewritten in a way 
that /i = A = 1, while in the quantum theory, this is not 
possible. 

We compared some of the results from these kink- 
antikink collisions with the heavy ion collisions such as 
are currently carried out at the RHIC in Brookhaven 
and in the future at the LHC at CERN in Geneve. Al- 
though the scalar i/)^ theory is too simple to serve as a 
toy model for this, it is interesting to see that we can cer- 
tainly compare some of the ingredients. The kinks then 
describe the colliding hadrons or ions, while the quasipar- 
ticles describe the formed "plasma". After the collision. 



most of the energy remains in the receding kink-antikink, 
while the resulting field shows a central plateau, flat in 
space-time rapidity, indicating hydrodynamic scaling. Its 
height only weakly depends on the initial speed of the 
kinks and scales with the coupling in a very straightfor- 
ward way. 

We found that the central region first has a finite speed 
of sound, of about 0.8 for A/m^ = 1/1.25 and about 0.4 
for X/m^ = 1/12. After a short while this goes over in 
a free expansion, with vanishing pressure and speed of 
sound. 

In this aspect the qualitative difference between the 
Hartree and classical approximations is small. Although 
it is difficult to compare the two at strong coupling at the 
later stage, as the Hartree approximation is "haunted" by 
the "symmetric minimum" at (/) = 0, at the early stage a 
meaningful comparison can be made. It shows very simi- 
lar results between the Hartree and classical approxima- 
tions. At smaller coupling the two are very similar even 
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during the later stage, the Hartree result shows somewhat 
more interaction at times around 50 — 100 and related to 
that a somewhat smaller energy density. 

It would be interesting to have some handle on the er- 
ror induced by the first order phase transition predicted 
by the Hartree approximation. Although the lack of a 
systematic expansion parameter in this approximation 
makes it difhcult to make a quantitative estimate, the 
above results seem to indicate a qualitatively small influ- 
ence, even at strong coupling, as long as the field is not 
moving too close to the artificial "symmetric minimum" . 
However, at the later stage in the strong coupling simula- 
tion of Fig. the field is actually fluctuating around this 
artificial minimum at (/> = and the results can no longer 
be trusted. Similarly, it is not expected that the Hartree 
approximation can accurately describe the dynamics of 
going through the phase transition, see Ref. 0| ■ 

In Ref. [SQl the Hartree "plasma" was also found to be- 
have as a "relativistic plasma," with speed of sound close 
to 1, similar to what we found for tm ^ 8. However the 
result in [s^l was obtained from a disintegrating Gaussian 
wave packet in the "symmetric ■rfiase," in which the ef- 
fective coupling is much smaller £|, instead of a colliding 
kink-antikink in the "broken phase." No free expanding 
phase was found in that reference. 



Finally, we have briefly looked at the connection be- 
tween a thermal Bose-Einstein distribution and the cre- 
ation and annihilation of kinks in the system. We found 
that one should only consider the complete field in the 
description of kinks, not just the mean field. This means 
that we have to look at quantum and ensemble averages 
only; the separate realisations give some impression of 
what is happening but do not describe the full theory. It 
is encouraging to see that although the kinks disappear 
from the mean field, our rough kink indicator does not go 
to zero, but becomes constant, i.e., the creation and an- 
nihilation rates become equal. However, we do find fewer 
pairs in the Hartree approximation than in the classical 
theory. This might be related to the higher instability of 
Hartree kinks, which can radiate quantum particles, in 
the Hartree description, more energy is carried by radia- 
tion than in the classical theory. 
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